Sum rule analysis of Umklapp processes and Coulomb energy: 
application to cuprate superconductivity 
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The third moment frequency sum rule for the density-density correlation function is rederived 
in the presence of Umklapp processes. Upper and lower bounds on the electron-electron Coulomb 
energy are derived in two-dimensional and three-dimensional media, and the Umklapp processes 
are shown to be crucial in determining the spectrum of the density fluctuations (especially for 
the two-dimensional systems). This and other standard sum rules can be used in conjunction with 
experimental spectroscopies (electron-energy loss spectroscopy, optical ellipsometry, etc.) to analyse 
changes of the electron-electron Coulomb energy at the superconducting transition in cuprates. 



I. INTRODUCTION. 



Theoretical progress in the understanding of micro- 
scopic origin of high-temperature superconductivity ap- 
pears to be ambiguous. The novelty and difficulty is to 
describe the strong effects of the electron-electron inter- 
actions which determine the strongly correlated phases 
of cuprates at various dopings: the Mott antiferromag- 
netic insulator, the "anomalous" metallic state and the 
superconductor. In fact, the fundamentally new micro- 
scopic origin of superconductivity(SC) in these materials 
is, perhaps, due to electron-electron interactions unlike 
standard phonon- mediated superconductivity We 
explore the general aspects of such a scenario with the 
help of sum rules for the density-density correlation func- 
tion. 

The particular question of interest is the origin of the 
condensation energy, the difference in the energy between 
the "normal" state extrapolated ioT ~Q K and the su- 
perconducting ground state. The famous BCS theory of 
superconductivity ^ is based on the assumption that 
the attractive interaction between electrons arises from 
the lattice vibrations. Here we would like to investigate 
a general alternative to the BCS phonon-mediated su- 
perconductivity, namely, that the superconducting state 
(e.g. the condensation energy) is promoted either by the 
long-range part of electron-electron Coulomb interaction 
or by the short-range part of electron-electron and static 
electron-ion interactions. The identification of the part of 
the full electron-ion Hamiltonian responsible for the con- 
densation energy would be an important step towards a 
complete and consistent theory of high-temperature su- 
perconductivity. 

Many proposals for the condensation energy have been 
suggested [^,||, most of them (except Ref. |3|) being 
based on phenomenological Hamiltonians. Experimen- 
tal confirmation of the origin of the condensation energy 
based on a phenomenological Hamiltonian would still 
require the validation of the phenomenological Hamil- 



tonian (by various other experiments) in order to con- 
struct a self-consistent theory. Alternatively, the essen- 
tial terms of the full |5j original electron-ion Hamiltonian 
(and changes of the expectation values thereof upon tran- 
sitions) can be established consistently from experiments 
with the use of sum rules at the outset. Such an ap- 
proach, taken in this paper, can potentially identify the 
origin of the condensation energy and the phenomeno- 
logical Hamiltonian sufficiently to describe superconduc- 
tivity and other strongly correlated phases. 

If the ion kinetic and ion-ion Coulomb energies are as- 
sumed irrelevant (or in other words, these terms do not 
change upon the phase transitions and do not determine 
the important correlated phases), the full [|| electron- ion 
Hamiltonian can be reduced to the following form 
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is the total den- 



where pg = J2k,a 4-q,aCk,a = Er, 6* 

sity operator. N is the number of electrons, and Q is 
a total volume. The first and second terms are the ki- 
netic and Coulomb interaction energy of the electrons. 
The third term describes the interaction of the electrons 
with the periodic potential of the lattice, which can be 
represented by the Umklapp pseudopotential with 
the sum over corresponding wavevectors k of the recipro- 
cal lattice 1^. The term of interaction between electrons 
and positive homogeneous ion background is omitted. In 
spite of making the "static lattice" assumption (the dy- 
namic lattice effects (e.g. phonons) are neglected), the 
Hamiltonian (|l|) is quite general. For instance, the Hub- 
bard model (and multi-band versions of it) is only a re- 
duced version of the Hamiltonian (^ which neglects the 
long-range part of the Coulomb interaction. Presumably, 
the Hamiltonian (|l|) is sufficient not only to describe the 
Mott insulating state and the metallic state at high densi- 
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ties (far away from the half-filled band) but also all other 
important phases of the cuprates. 

We can analyse the electron Coulomb energy in situa- 
tions of different dimensionality. In an isotropic medium 
the three-dimensional (3D) Coulomb potential is Vq = 

2 

where eoo is the high-frequency dielectric con- 
stant due to the screening by the core electrons In 
a two-dimensional (2D) plane, the Coulomb potential 



is Vn 



For a layered electron gas with inter- 
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plane distance d (relevant for the discussion of single- 
layer cuprates) the Coulomb interaction is 



^ e^d sinh(qd) 

2eoeoo'Z cosh{qd) — cos{qzd) 



(2) 



with different dependencies on the wavenumbers q paral- 
lel and qz perpendicular to the planes The Coulomb 
potential of a layered gas becomes three-dimensional 

Va 



^ in the long- wavelength limit (qzd <ti 1 and 

9 

for short 



qd ^ I) and two-dimensional, Vq 



wavelengths {qd 3> 1, independently of q^ momentum) 



II. SUM RULES. 

The electron-electron Coulomb energy can be related 
to the density response function x(g,w), and thus it is 
instructive to analyse various sum rules for this response 
function In particular, the expectation value <Vc> 
of the Coulomb energy can be written in the form 



< Vc,, >= ^Vq J —Imx{q,u^)coth ^— j , 



(3) 



where n = N/fl is the density of the electron system, 
and the density-density correlation function x((j, w) pC| ] 
is defined in the standard way [|liyi2[ 



■ p+co 

X{q, dte<-+^'^' < [p{q, t),p+{q, 0)] > . (4) 

< Vc.q > is the expectation value of the partial Coulomb 
energy corresponding to a particular value of momentum 
q. The structure factor Sq 



So 
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< PqP-q > -NSq^o, 



(5) 



can be expressed again through the imaginary part 
Imx{q,uj) using the fluctuation-dissipation theorem 
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Various sum rules for the imaginary part of the suscep- 
tibility Imxiq, '-^) (valid for arbitrary dimensional system 
with appropriate form of Coulomb potential Vq) can be 
derived by calculating the commutators of the density op- 
erator and the Hamiltonian. To calculate the well-known 
/-sum rule (or the first moment sum rule) it is neces- 
sary to calculate the expectation value of the operator 
[[pq, H], p-q\. Another additional sum rule is the causal- 
ity (Kramers-Kronig) relation. These two well-known 
sum rules are |12l 
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To calculate the w-^-moment sum rule the expecta- 
tion value of the operator [[[[p^, i/], /5_q] must be 
calculated. The result of a long calculation is 
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is the 



projected kinetic energy operator, and q ~ q/\q\ is a unit 
vector along the direction of q. This higher order sum 
rule (Eqn. |9|) is convenient for analysis of high frequency 
transitions in the density response, because it weights 
higher frequencies by a factor w^. The existence (or con- 
vergence) of the third moment sum rule can be demon- 
strated by showing that all terms on the right-hand side 
of Eqn. ^ are finite. The first and the third terms are 
expected to be finite, because they are essentially related 
to the finite expectation values of the electron-lattice and 
the kinetic energy in the ground state. The convergence 
and g-dependence of the last term will be discussed in 
detail separately for 2D and 3D cases. 

The three sum rules (Eqn. 0-|^) allow us to derive upper 
and lower bounds p5[ on the electron Coulomb energy, 
and hence to discuss the possible changes of Coulomb 
energy due to the phase transition. Since the imaginary 
part of the susceptibility Imxiq, ^) a real positive def- 
inite function, the two Cauchy-Schwartz inequalities can 
be written for the partial Coulomb energy < Vc >q at 
T = K 
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It is convenient to introduce the notional "plasma fre- 

, 2y \ 1/2 

quency" ujp{q) defined by ujp{q) — I ) . In fact, 

defined as above the "plasma frequency" has the right 
asymptotics for the corresponding plasma waves in the 
three-dimensional case (ujp ^jj{q) — e^n^D / {^o'm)) and 
the two-dimensional case (1^^20(9) ~ e^n2D'?/(2eom)). 
Another relation, which is useful in order to rewrite the 
J_i sum rule for (7 — s- 0, can be derived in the long- 
wavelength limit [q <C qxF, where qxF is the inverse of 
Thomas-Fermi screening length) for the full susceptibil- 
ity x(g, 0), if we express x(g, w) through the "bare" (or 
local) susceptibility xo{(1t^)'- 



X(9,^) 



Xo('Z,t^) 



1 + VqXoiq,^^) 



(11) 



(this is an exact result, not an RPA approximation, pro- 
vided xo is defined in terms of the relevant irreducible 
diagrams). The "bare" susceptibility xo(<Z, 0) of the elec- 
tron liquid is assumed to be finite (see section 4.1 of Ref. 
|l2[ for the compressibility sum rule) , then for g <C (?tf 



x{q,^ = 0) = 
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The discussion up to this point is valid for any dimen- 
sionality of the system (with the corresponding form of 
Coulomb potential Vq). In what follows we analyse the 
two-dimensional and three-dimensional cases separately 
and find important differences. 

We consider first the 3D case. The leading terms at 
small q for the third moment sum rule are 



^ < ^ > +r— 



(13) 



where A = ^Y^jK^)^i~l)U-^p^ |T|]. The third and 
fourth terms in Eqn. |^, being proportional to q'^ and q^ 
powers, are subdominant. The last term in (pf) 
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{pq + q'Y - Vp{pq)']nSp 



(14) 



requires careful analysis of its convergence and q- 
dependence. The g-dependence (and convergence), due 
to the part of the integral over small momenta p, is evi- 
dent, since the upper bound on the pair-correlation func- 
tion Sp (see Eqn. ||, ^ is Sp < Bp^^+^'>/^ (where B is a 
constant, D is dimensionality) 



dSVpipq)^[Sp+q -Sp]^ / dSVp{pq)^Sq (15) 



for small momenta p < q. Since in 3D case the upper 
bound Sq < Bq^, the contribution from small momenta 



p integration is at least of order q'^ (or higher power of 
q for small q). To analyse the Q'-dependence of the part 
of the integral (Eqn. |lj) from the integration over large 
momenta p, we use the "cusp theorem" JT8| , which gives 
the asymptotic behaviour at large momenta p (3D case) 



Sp = l — + o{ — ), 

p-i pO 



(16) 



where C3D is some constant. The integral of Eqn. [ij with 
5p = 1 is identically zero (this is why Sp can be substi- 
tuted by S'p — 1 if convenient). It is easy to see that the 
integral over p with the second term of expansion C^u jp^ 
is convergent. By expanding in powers of (g/p) for large 
p, the leading g-dependence of the considered term 



d^pVpm\Sp^, -Sp]= I rf3 /p(Pyg3D ^ 



(1 + [q/pf + 2{q/p)cos6Y 



(17) 



is found to be proportional to q'^. The terms, propor- 
tional to and other odd powers of q, are equal zero 
after the integration over angle 9 {cosO — {pq)/pq is the 
angle between p and q)). It can be seen that for large 
q the last term (Eqn. |lj) grows no faster than q^ as 
well. What is important for the ensuing discussion is 
that, both for small and large g, the last term of Eqn. ^ 
has subleading 5'* dependence on the wavevector q. 

The upper and lower limits on the partial Coulomb 
energy can be conveniently written in terms of the 3D 
"plasma frequency" 



+ oiq'). 



(18) 
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The inequalities ( p^ ) have an interesting significance for 
the possibility of gaining energy, in the SC condensation, 
from small-g modes. We first notice that in 3D case if 
< ^ >= in both normal and superconducting states, 
which is certainly the case if there is no crystalline po- 
tential, the terms proportional to q'^ (the third and last 
terms of the right-hand side of Eqn. ^ determine the 
difference between the upper and lower bounds. There- 
fore if < A >~ 0, then for given q the maximum pos- 
sible saving is proportional to {h/2)LUp,3jj{q/qQ)'^ (where 

(Jo ~ qTF,3D, qTF,3D \[^^^ in 3D, kp is the Fermi 
wavevector), which is a negligible portion of the partial 
Coulomb energy ^ujp^so at long wavelengths. In fact, 
in the absence of the Umklapp processes the sum rules 
(Eqn. H and ^) essentially fix the density spectrum at 
long wavelengths q < qrp io the plasma pole contribu- 
tion: Imx{q,Lo) ~ ^ft"p^o(g) _ which sat- 
isfies completely all three sum rules at g ^ . Other 
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terms in the J3 sum rule become comparable with the 
dominant term q only at q > qrp- Therefore, a 

large saving of Coulomb energy in 3D case in the ab- 
sence of the lattice due to the phase transition is possible 
only for short wavelengths q > qTF,3D- 

In the presence of the periodic lattice potential 
(< A 0), it can be seen from the inequalities ( p^ that 
for q < qTF,3D the maximum theoretical saving is a fi- 
nite fraction of the zero-point plasma energy ^ujp^3D{q). 
Therefore, substantial saving of the Coulomb energy in 
the long-wavelength limit is possible only in the pres- 
ence of strong crystalline potential. In 3D (where the 
Coulomb interaction Vq = e^/eog^), the Umklapp term 
< A > contributes in the same leading order of pow- 

ers of q into the J3 sum rule as the Coulomb term q ^Vq 
(see Eqn. The Umklapp term is then responsible for 
the finite width (or "lifetime") of the plasmon peak. 
In the 2D case, the leading terms are again 



J3 = ^ <A> +q^^Vq^2D- 



(19) 



Other terms are subleading and proportional to q^ and 
q^ powers. The last term can be analysed similarly to 3D 
case and be shown proportional to q* at small q. Using 
the "cusp theorem" for 2D case [fsl 



Sn 
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+ — ), 



(20) 



we can expand in powers of q/p for large momenta p and 
keep the leading term as a function of q/p: 



J d?pVp^2D {pq)^[Sp+q - Sp]~ J (fp 



(1 + (g/p)2 + 2{q/p)cos0)y^ 



2 Vp^2D{pq}'^C2D 



Due to the angle integration, terms proportional to q'^ 
and other odd powers vanish. For large q, the last term 
grows no faster than also (due to the first term in the 
bracket of Eqn. |l^. In the absence of Umklapp scat- 
tering (< A >— 0), the density spectrum is given by 
the expression Imx{q,uj) ~ '^^^y^"^'^^ 5{uj — ujp^2D{q)), 
which satisfies all three sum rules at small q. The max- 
imum possible saving of partial Coulomb energy at long 
wavelengths is of order |t^p,2D(g) q^l^o (^^'^re qTF,2D — 

2Tr r^h'^ ^ 2^ Thomas- Fermi screening wavevector). 

2 

The presence of the Umklapp term ^ < A > has a 
much more dramatic effect on the density spectrum in 
2D, because this term, proportional to q^, has a lead- 
ing power of q in the third moment sum rule at small 
q dominating over the Coulomb term (unlike 3D case, 
where the Umklapp term has the same power-g depen- 
dence as the Coulomb term). The density spectrum can- 
not be even approximated by a plasma pole expression. 



and the plasmon is never really a well-defined excitation 
in 2D in the presence of Umklapp scattering. It means 
that the Umklapp scattering modifies strongly (or "non- 
perturbatively" ) the spectrum of the density fluctuations, 
and the spectrum is dominated by multi-pair and pair 
excitations rather by plasmon. Of course, the mere ex- 
istence of a large value of < 4 > is in itself perfectly 
compatible with a traditional textbook picture, in which 
the sum rules are satisfied by taking proper account of 
interband transitions; in such a case there is no a pri- 
ori reason why a plasmon associated with the excitations 
of the conduction band must automatically be ill-defined. 
However, in a system where Umklapp (quasi-momentum- 
nonconserving) interactions between the Bloch quasipar- 
ticles are strong it seems natural that these alone could 
give rise to a substantial value of Imx even below the 
first band gap; particularly in view of the above remark 
about the enhanced effect in 2D, it is tempting to view 
the so-called midinfrared peak in the cuprates in this 
light. 

The upper bound on the partial Coulomb energy is still 
the half of the plasmon energy {ujp^2D{q) ^ y/q), but the 
lower bound at small g if < A >7^ is essentially given 

by 



< Vc.q > 



> 



^p.2D{q) 



(22) 



<A> 



and so in the limit g — > 0, 



"^,20(9) 



1/2 



<A> 



q. (23) 



Therefore, bounds, based on sum rules, are compatible 
(21) with saving of almost all Coulomb energy -1^^,20 (9) in 



2D when < A>y^O. 

It is necessary to mention extensive literature (for in- 
stance, Ref. |l^ ) using sum rules (in particular, third 
moment sum rule) in order to analyze and derive various 
local-field corrections and approximations of the density 
response, whereas our goal in this paper is to analyze 
general constraints on the electron Coulomb energy at 
small q without relying on any approximation. 

It is also interesting to discuss briefly for comparison 
the interaction energy of a many-particle system inter- 
acting via a short-range potential {Vq — > Vq = const, 
for g — !■ 0). The upper bound on the partial interac- 
tion energy < Vint,q > is given by the "acoustic mode" 

|(^o^i^xo(g,o))V2 



Von 
ms 



q > < Vnt.,q >, 



(24) 



where Xo('?jO) is given by the compressibility sum rule 
]l2[: Xo('Z, 0) = where s is the velocity of 
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sound. Therefore, the maximum available interaction en- 
ergy at long wavelength is insignificant (especially when 
weighted by the phase volume). For instance, it implies 
that in most phase transitions in neutral systems (i.e. 
many-particle systems interacting via a short-range po- 
tential) the interaction energy is saved predominantly at 
short distances. 



III. EXPERIMENTAL PROBES. 

The remainder of this paper is devoted to a brief dis- 
cussion of the experimental spectroscopies which should, 
at least in principle, be able to shed light on the origin 
of the condensation energy in the SC transition, and the 
inferences which we may currently draw from them (cf. 
also Ref. section 4.2). For simplicity we will consider 
explicitly a single-plane cuprate such as T/ — 2201, so that 
in the normal phase there is only one characteristic length 
(other than, possibly, the electron mean free path) large 
compared to the quantity q^p, namely the interplane 
spacing d (typically ~ 10 A, i.e. ~ 10— 20(7y^)(note nota- 
tional differences from Ref. |^). The case of a multi-layer 
cuprate such as Bi — 2212 is more complicated, as there 
is now a second "large" characteristic length, namely the 
intra-bilayer spacing ('^ 3 — 5 A); however, the general 
pattern of the results is unchanged. In addition, we will 
assume tetragonal symmetry. 

The two spectroscopies which most directly probe the 
Coulomb energy, or something closely related to it, are 
electron energy loss spectroscopy (EELS) ||l^ and optical 
reflectivity |^^; in the latter case we shall assume that 
ellipsometric measurements are possible in the interest- 
ing frequency regime |21 23] so that we may deduce the 



relevant complex dielectric constant without the use of 
Kramers- Kronig relations. In a bulk isotropic 3D metal 
the situation is very simple: to the approximation that 
we neglect multiple scattering and the effect of the ionic 
cores, the transmission EELS cross-section a{q, lu) is a di- 
rect measure of the quantity V^Imx{q, i^), where x{Qj ^) 
is the "true" density susceptibility as defined in Eqn. ^. 
Since in 3D the longitudinal dielectric constant e||(q,(jj) 

is identically equal to 1 + -^^Xoily'-^) with the help of 
Eqn. ^l], we find the simple result 



a(q,uj) = const^Im 



=11 (9, w) 



(25) 



where the constant is of purely geometrical origin and can 

is usually 



11(9,'^) 



be calculated, and the quantity Im 
known as the loss function. This formula is valid for arbi- 
trary q, including values larger the inverse lattice spacing 
in the approximation of neglect of the direct scattering 
effect of the ionic cores. If the latter is taken into account, 
the effect is to multiply the formula ( p5| ) by a factor which 
is in general a function of q but not of oj provided that 



the latter is small compared to typical core excitation 
energies (~ 20eV), and expected to be unaffected by the 
superconducting transition. Writing out the integrand of 
Eqn. ^ explicitly in terms of xoilj^) (see Eqn. [Tl| ) and 
using the 3D relation between the latter and e\^{q,uj), we 
see that apart from a function of q the transmission EELS 
cross-section is a direct measure of the Coulomb energy 
locked up in dqdio. 

It is well known that the condensation energy due to 
the SC transition is extremely small (of order 10"** eV 
per electron (or per unit cell)) in comparison with the 
atomic energies (10 eV). It implies stringent require- 
ments on experimental techniques, nevertheless, changes 
associated with the SC transition were observed by optics 
at mid-infrared frequencies ||2^,|2^ of magnitude sufficient 
to provide the condensation energy (measured directly 
by specific heat measurements). It is worth noting that 
we do not discuss the changes at frequencies compara- 
ble or lower than a superconducting gap (although these 
changes of course are most remarkable consequences of 
superconductivity!), because the change of Coulomb en- 
ergy associated with this region of frequencies is negligi- 
ble (if limited to small momenta q <C qxp)- The optical 
refiectivity measurements 21,2^] have enough precision 
to explore the type of questions discussed in this paper, 
while it is hoped that the transmission EELS can achieve 
required accuracy in the near future. 

In a 3D bulk metal ellipsometric optical measurements 
can measure the complete transverse dielectric constant 
e±{q,LLj), in the limit (f — > 0, and hence the correspond- 
ing "transverse" loss function Im r — 7 . Since in 

o [ ej_(i3,w)J 

the normal phase, at least, there should be no distinc- 
tion, in the limit q — > 0, between ej[((7, w) and e±{q,uj), it 
follows that in this phase the (f — > limit of "loss func- 
tions" measured by EELS and by optics should coincide. 
It is a somewhat delicate question, once one renounces 
reliance on some specific model such as the Fermi liquid 
one, what is to count as "the g — > limit" ; in addition to 
the obvious scale qrp or qp, it is not immediately clear 
that the inverse electron mean free path l/l might not 
be a relevant quantity. However, it is plausible that this 
quantity should not play a major role for uj in the mid- 
infrared region, so that we shall tentatively take "q — > 0" 
to mean in the 3D case q <ti qxF, qp- 

Some care is needed in adapting the above results to 
the case of a layered material such as cuprates, even if 
we specialize (as we shall) to the limit q^ ^ l/d, which 
is automatically fulfilled in optical experiments and may 
be satisfied in transmission EELS by a suitable choice 
of geometry. In the context of EELS experiments, we 
now have to distinguish the cases qd ^ 1 and 3> 1 
(where q is the ab-plane component of the momentum 
loss). In the former case the 3D bulk formula (Eqn. [25| ) 
applies unchanged, provided that e\\ {q, uj) is defined to be 
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the tensor component of the longitudinal dielectric con- 
stant corresponding to current flow in the ab-plane; note 
in particular that, at least in the normal phase, we ex- 
pect that Im 
the limit q — > 



0. 



is nearly independent of Iff] in 

In the opposite limit qd ^ 1, we could 
choose to continue to use the 3D formula(Eqn. ^5[), but 
would than find that the e|| {q, uj) so defined has a strong 
explicit dependence on q. A much more natural con- 
vention in this limit is to treat the scattering as occuring 
independently from the different Cu02 planes, and to de- 
fine a two-dimensional ( "per-plane" ) version x^\qt^) of 
XoC'Ziw), or equivalently a quantity (cf. Ref. Q), section 
4.1) 



K{q,L0) = 



1 



xl3\<l^^)^^iHi<l,^)-^b), (26) 



2eor 



where 611(5,0;) is the "natural" definition of the 3D bulk 
ab-plane dielectric constant, i.e. the quantity which re- 
lates the local polarization to the local field, and is ex- 
pected to be nearly constant over a range q 3> l/d, 
and eb is the "background" (off-plane) contribution to 
it (cf. Ref. 1^). With this definition we find that, apart 
from factors depending only on q, both the transmission 
EELS cross-section and the (single-plane) Coulomb en- 
ergy locked up in the range dqdco are proportional to the 
quantity 



Im 



1 



1 + qK{q,u))/es 



(27) 



where esc is the dielectric constant which screens the 
Coulomb interactions of the in-plane electrons (note that 
in general esc is not equal to e^). Thus, just as in the 
bulk 3D case, the transmission EELS cross-section is a 
direct measure of the Coulomb energy locked up in the 
relevant region of {q, ti;)-space; note however that Eqn. |2^ 
introduces an extra explicit g-dependence which absent 
in the bulk case. This subtlety seems to have been over- 
looked in the analysis of existing normal state EELS data 
[ pof on the cuprates, where it seems to be assumed that 
even in the regime qd ^ 1 EELS experiments measure 
the "bulk" e||(g,w). In the intermediate case {qd ^ 1) a 
similar analysis using Eqn. ^ is possible, but will not be 
given here. 

In the case of optical experiments on the cuprates 
(with the sample surface assumed to lie in the ab-plane) 
we always have <C 1, ^ 1 , and thus at first 
sight we would expect the complex dielectric constant 
inferred from reflectivity measurements to be identical 
to the e||(q,cj) inferred, via Eqn. from EELS ex- 

periments in the limit q —^ 0. Existing measurements in 
the normal phase do appear to be consistent with this 
prediction. In the superconducting state, however, there 
are three complications: first, it is not clear that even in 
the mid-infrared regime the Cooper-pair radius ^0 is not a 



relevant length scale, so that it may be illegitimate to use 
the "true" q^O {q^o ^ 0) behaviour observed in the op- 
tics to infer the behaviour in the regime l/^o q ^ qxF- 
Secondly, it is not completely obvious, particularly in the 
former limit q^o ^ 1, that the finite- frequency longitudi- 
nal and transverse dielectric constants must be equal in 
the superconducting state. This latter complication may 
be somewhat mitigated by a third consideration, namely 
that at the non-normal angles of incidence necessary in 
the ellipsometric technique what is measured, in a layered 
material, is not simply e± but a combination of e± and e|| ! 
We will not attempt to develop those points further here, 
but will rather use them to draw the conclusion that, 
while the spectacular changes observed 22| in the op- 



tically measured dielectric constants of the cuprates at 
and below the superconducting transition are strongly 
suggestive, a quantitative test of any scenario (such as 
the mid-infrared one of Ref. 1^) which attributes the en- 
ergy saving largely to a regime of q small compared to 
qTF but large compared to ^0 (and 1/d) will require accu- 
rate transmission EELS data taken across the transition, 
something which (as regards the mid-infrared regime of 
frequencies) does not to our knowledge at present exist. 

We finally address head-on the question: which of the 
three terms in Eqn. |l| is (are) reduced in the supercon- 
ducting transition? If it is the second (Coulomb) one, in 
what regime(s) of q and cu does the saving predominantly 
occur? Part of the interest of this question is that as we 
have seen above, a conjectured answer can be tested di- 
rectly in transmission EELS experiments. 

We start by recalling a well-known result: since the 
original ("true") Hamiltonian of the N-body system 
(the nonrelativistic limit of the Dirac Hamiltonian) is 
composed exclusively of kinetic energy terms and (un- 
screened) Coulomb interactions, the virial theorem im- 
mediately tells us that the change in total kinetic energy 
(of electrons and ions) must be exactly minus half that of 
the total Coulomb energy (electron-electron, electron-ion 
and ion- ion) , and thus Coulomb energy must be saved in 
the superconducting transition (and indeed in any other 
phase transition into a lower-energy state). While this 
conclusion is very generic and rigorous, it is not usually 
regarded as giving much insight into the "mechanism" 
of superconductivity in the cuprates (or for that mat- 
ter in the classic superconductors) since the term "mech- 
anism" is often held to refer to a low-energy effective 
Hamiltonian in which the separation of the original ki- 
netic and potential energies may no longer be explicit. 
However, "intermediate-level" effective Hamiltonian (|l|) 
is sufficiently close to the original truly first-principles 
one that the virial-theorem result for the latter might at 
least suggest that it is one or both of the last two terms 
of (|l|) which are saved |2^ . 

If we assume for the sake of argument that Coulomb 
energy is indeed saved, then where in the space of q and 
oj is it saved? It is at this point that the sum rules argu- 
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ments of section 2 come into their own. For convenience 
we reproduce here the three relevant sum rules with terms 
of relative order q"^ and higher omitted on the right-hand 
sides: with notation as in Eqns. 0,||J|; 

J-i^^. (28) 



Ji = 



m ' 

9 4 2 

n <l \, q i 
-^Vq + < A > 



(29) 
(30) 



where in the case of a layered system Vq is given by 
Eqn. H, and tends to e^/(eoeooQ^) for <C 1 (and 
<C 1, see above) and to e^/(2eoeoo<z) for » 1, 
and where A is defined below Eqn. The {T = K) 
contribution < Vc^q > to the expectation value of the 
Coulomb energy from wave vector q is, up to a factor, 
just Jq (see Eqn. ||). 

The arguments of section 2 now show that the maxi- 
mum change of < Vc^q > in the "essentially 3D" regime 
qd ^ 1 is proportional to a finite fraction of the "3D plas- 
mon energy" but weighted by q^ of the phase space (d^q), 
and hence the maximum total saving possible from this 
regime is very small. On the other hand, the contribution 
from the regime 1/d <^ q ^ qxF (where the truncated 
forms Eqns. 28-30| are still a good approximation) can 
be of order q^ (cf. the conclusion after Eqn. ^), pro- 
vided the quantity < A > is substantial, while the phase 
space allows significant saving of Coulomb energy since 
dq^d^q ~ {2Tr / d)2Trqdq and 1/d ^ q ^ qrF- Thus, in a 
quasi-2D system with a large value of < A > substantial 
energy is in principle available for saving in this small-g 
regime. To estimate the value of < A > we return for 
a moment to the limit qd ^ 1 and refer to the normal 
state optical loss function data: using Eqns. p8[]30| with 
the appropriate form (e^/(eoeoo?^)) of Vq, we see that 
the quantity < A >„ is given, in the natural units of 
n^e^/(eoeoo) by the expression ("n" is the normal-state 
value) 



<A>n 



= .h - Jl/J- 



(31) 



Although a strict evaluation of the right-hand side of 
Eqn. ^from the optical loss- function requires us to know 
about the effective frequency cutoff (since at high fre- 
quencies there will be contributions to e{q, lj) from "core" 
processes not described by Eqn. it is clear that the 
mere existence of mid-infrared(MIR) peak extending over 
an order of magnitude in frequency already implies that 
it at least of order of 1. Thus, a very appreciable fraction 
of the Coulomb energy locked up, in the normal state, 
in the \ow-q, MIR-frequency regime is in principle avail- 
able for saving in the SC transition(or indeed in other 
possible phase transitions). Whether it is in fact saved 
and to what extent, as in fact postulated in the "MIR 



scenario" of Ref. 1^ , depends of course on the cost of the 
formation of the Cooper pairs in kinetic and/or static 
lattice energy. Actually, rather than asking as above for 
the fraction of the Coulomb energy which is in principle 
available for saving (something which is not that signifi- 
cant if the original value is itself small), it may be more 
informative to estimate the relative contribution of the 
small-g regime in 2D and 3D for a given change in xo 
due to the phase transition. Taking into account both 
the phase space factor and the extra factor of q in the 
denominator of the expression ( p7[ ) in the 2D case, we 
find that in the regime where the {qK) dominates the 
contribution of small q is proportional to q^ in 3D but 
to a constant in 2D, so that the relative importance of 
the long-wavelength regime is vastly enhanced in the 2D 
case. 

On the experimental front, it has to be said that as 
noted in Ref. |2^, the optical data, if extrapolated into 
the relevant (g^o ^ 1) regime with several other as- 
sumptions, indicate rather the opposite, i.e. that the 
Coulomb energy associated with the MIR regime actually 
increases in the SC state. However, because of the var- 
ious considerations noted above, this extrapolation may 
be problematic, and a definitive test of the MIR hypothe- 
sis must await quantitative transmission EELS measure- 
ments across the superconducting transition (or a better 
theoretical understanding of the generic q-dependence of 
e{q,U)) in the SC state). 

In sum, we analyzed the electron-electron Coulomb en- 
ergy in the presence of the periodic lattice potential us- 
ing various sum rules for the density-density response 
function. We believe that in this paper we have made it 
plausible that two specific properties of cuprates, namely, 
(a) the layered (two-dimensional) structure of the Cu02 
planes and (b) the occurrence of a broad and strong peak 
in the optical loss function, can be essential ingredients 
in the occurrence of high-temperature superconductivity 
in these materials by conspiring to the saving of small-g 
{q < qrp) part of electron-electron Coulomb energy. 
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and Neil Ashcroft. A.J.L. acknowledges the support of 
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